CT invariant quantum spin Hall effect in a ferromagnetic graphene 
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We predict a quantum spin Hall effect (QSHE) in the ferromagnetic graphene under a magnetic 
field. Unlike the previous QSHE, this QSHE appears in the absence of any spin-orbit interaction, 
thus, arrived from a different physical origin. The previous QSHE is protected by the time-reversal 
(T) invariance. This new QSHE is protected by the CT invariance, where C is the charge conjugation 
operation. Due to this QSHE, the longitudinal resistance exhibits quantum plateaus. The plateau 
values are at 1/2, 1/6, 3/28, ... , (in the unit of h/e 2 ), depending on the filling factors of the spin-up 
and spin-down carriers. The spin Hall resistance is also investigated and is found to be robust 
against the disorder. 
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In the years since the spin Hall effect (SHE) has 
been discovered, it has generated great interest.— ~— In 
the SHE, an applied longitudinal charge current or volt- 
age bias induces a transverse spin current due to the 
spin-dependent scatterings 2,3 of the spin-orbit interac- 
tion (SOI)l Soon afterwards, the quantum SHE (QSHE) 
is also predicted^ The QSHE occurs in a topological 
insulator in which the bulk material is an insulator with 
two helical edge states carrying the current^ The edge 
states, with opposite spins on a given edge or opposite 
edges for a given spin direction containing opposite prop- 
agation directions, lead to a quantized spin Hall conduc- 
tance. The QSHE is a new quantum state of matter 
with a non-trivial topological property. The existence 
of QSHE was first proposed in a graphene film in which 
the SOI opened a band gap and established the edge 
states.— £ But the sequent work found that the SOI in 
the graphene was quite weak and the gap-opening was 
small, so the QSHE was difficult to observed Soon af- 
terwards, the QSHE was also predicted to exist in some 
other systemsjiS^— Recently, the QSHE was successfully 
realized in the CdTe/HgTe/CdTe quantum wells, and a 
quantized longitudinal resistance plateau was experimen- 
tally observed due to the QSHEii 2 . 

Another subject that has also been extensively inves- 
tigated in recent years is the graphene, a single-layer 
hexagonal lattice of carbon atoms^ after it has been 
successfully fabricated The graphene has a unique 
band structure with a linear dispersion near the Fermi 
surface, giving it many peculiar properties. For example, 
the quasi-particles obey the Dirac-like equation and have 
rclativistic-like behaviors, and its Hall plateaus are at the 
half-integer values. 

In this Letter, we predict a new kind of QSHE in 
a ferromagnetic graphene. Let us first imagine a two- 
dimensional system consisting of the following character- 
istics: (i) its carriers contain electrons and holes; (ii) both 
electrons and holes are completely spin-polarized with 
opposite spin polarizations. When a high perpendicular 
magnetic field is applied to the system, the edge states 
are formed and the carriers move only along the edges. In 



particular, the electrons (with their spins up) and holes 
(with their spins down) move in opposite directions on 
a given edge (see the inset in Fig. la). Therefore, the 
QSHE automatically exists in this system. Although the 
ordinary metals (or doped semiconductors) cannot meet 
the above two characteristics, a ferromagnetic graphene 
does. Recently, several approaches to realize a ferromag- 
netic graphene have been suggested J2r— For example, 
the ferromagnetic graphene can be realized by growing 
the graphene on a ferromagnetic insulator (e.g. EuO). 18 
For a ferromagnetic graphene, as soon as the Fermi en- 
ergy Ep is tuned to lie between the spin-up and spin- 
down Dirac points (see the inset in Fig. lb), the above 
two characteristics are met and the QSHE occurs. In the 
following calculations, we consider four- and six-terminal 
graphene Hall bars (see the insets in Fig. la). The results 
reveal that the transverse spin current and spin Hall re- 
sistance indeed show the quantized plateaus because of 
the QSHE. 

Comparing this new QSHE with the previously-studied 
QSHE, there are two essential differences: (i) The previ- 
ous QSHE comes from the SOI and the proposed systems 
all contain the time-reversal sy mmet ry^~— i 1 ' 1 1 while the 
present QSHE exists without the SOI and breaks the 
time-reversal symmetry. However, this new QSHE is pro- 
tected by CT invariance. (ii) In the previous QSHE, the 
edge states only carry a spin current while at equilibrium; 
in this QSHE system, the edge states carry both spin and 
charge currents at equilibrium with the two edges states 
being CT partners of each other (see the inset in Fig. la). 
Thus, this is a new kind of QSHE and the system is a 
new type of topological insulator. Due to the topological 
invariance, the plateaus of the spin Hall resistance are 
robust to disorder or impurity scattering. So the plateau 
is very stable and its value can be used as the standard 
value for the spin Hall resistance. 

In the tight-binding representation, the four- or six- 
terminal ferromagnetic graphene device (see the insets in 
Fig. la) can be described by the Hamiltonian: 20 
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FIG. 1: (Color online) The Hall conductance I e /V (a) and 
spin Hall conductance I s /V (b) vs. the energy eo for N = 
80 and <f> = 0.005. The two insets in (a) are the schematic 
diagram for the four- and six-terminal graphene's Hall bars. 
The inset in (b) is the schematic diagram for band structure 
of the ferromagnetic graphene while eo + M > Ef > eo — M. 



where and a ia are the annihilation and creation op- 
erators at the discrete site i. eo is the on-site energy 
(i.e. the Dirac-point energy), M is the ferromagnetic 
exchange split**, and t is the nearest neighbor hopping 
element. Here, the whole device, including the center re- 
gion and four or six terminals, is made of the ferromag- 
netic graphene. With the presence of a perpendicular 
magnetic field B, a phase factor (f>ij is added to the hop- 
ping element, <j>ij = J? A-dl/tfio with the vector potential 
A = (~By, 0, 0) and </> = h/e. 

The transmission coefficient T pqr7 (e) from the termi- 
nal q to the terminal p with spin a can be calcu- 
lated from the equation^. T pqa {e) = Tr[T pa G r a T qa G%}, 
where T pa (e) = i[E pcr (e) — £" (e)], the Green functions 
Gl(e) = [G«(e)]t = l/[e - H-« - £ p EJJ, and is 
the Hamiltonian of the center region. The retarded self- 
energy Y? va (e) due to the coupling to the terminal p can 
be calculated numerically^ After obtaining the trans- 
mission coefficient, the particle current in the terminal 
p with the spin er can be calculated from the Landauer- 
Bitttiker formula: I pa = (1/h) J de^2 q Tpq a (e)[f qcr (e) - 
f ptr (e)], where f p<7 (e) = l/{exp[(e - ii pa )/k B T] + 1} is 
the Fermi distribution function in the terminal p, with 
the spin-dependent chemical potential [i pa and the tem- 
perature T. In the following numerical calculations, we 
take t = 1 as the energy unit and only consider the zero 
temperature case (T = 0), as the thermal energy fc^T 
is normally much smaller than other energy scales in the 
problem. The sample width is denoted by N, and the 
insets of Fig. la show a system with N = 3. In the calcu- 



lations, we choose N = 80 and 40, and the correspond- 
ing widths are 33.9nm and I6.9nm. The magnetic field 
is described by the <f> with 4> = (3y / 3/4)a 2 i3/0o and the 
magnetic flux in a honeycomb lattice is 2cj). 

We first consider the four-terminal device (see the in- 
set at the top right corner of Fig. la) and a small bias V 
is applied between the longitudinal terminals 1 and 3 to 
study the induced charge current I ne [I ne = e{I n -\ + hi)] 
and spin current I ns [I ns = (^/2)(/„ t - 7 l4 )] in the 
transversal terminals 2 and 4. Here the boundary con- 
ditions for the four terminals are n\-\ — //ij, = eV/2, 
M2t = m = 0, (i3t = /i 3 j. = -eV/2, and = fi H = 0. 
The currents in the terminals 2 and 4 satisfy the rela- 
tions: I 2e = -he = h and I 2s = ~I is = I s . Fig. la 
and lb show the Hall conductance I e /V and spin Hall 
conductance I s /V versus the Dirac-point energy eo, re- 
spectively. For a non-ferromagnetic graphene (M — 0) 
under the high magnetic field (<f> = 0.005), I s /V is zero 
and I e /V exhibits the plateaus at odd integer values 
ne 2 /h (n = ±1, ±3, ...) due to the quantum Hall 
effect (QHE). These results have been observed in re- 
cent experiments 1 15 i 16 For a ferromagnetic graphene with 
M 7^ 0, the spin current emerges (see Fig. lb) since the 
QSHE. The spin Hall conductance I s /V also shows the 
quantized plateaus. By considering the edge state under 
the high magnetic field, the plateau values of I s /V and 
h/V can be analytically derived to be at (z/f — z^)e/87r 
and [y^ + v\)e 2 /2h^ where v a is the Landau filling fac- 
tor for spin a. In particular, when |erj| < |M|, in which 
case the Fermi energy Ep (Ep = 0) is located between 
the spin- up Dirac-point eo — M and the spin-down Dirac- 
point eo + M, I e is zero and a net quantum spin current 
emerges in the transversal terminals. In addition, if in 
the open circuit case, the spin accumulation emerges at 
the sample boundaries instead of the spin current^ 

Since the QSHE can give rise to quantum plateaus in 
resistances, we next study the longitudinal and Hall re- 
sistances in the six-terminal Hall device (see the inset 
in the lower left corner of Fig. la). Now we consider 
a small bias V applied to the longitudinal terminals 1 
and 4. The transversal terminals 2, 3, 5, and 6 are all 
voltage probes, their charge currents vanish (I pe = 0) 
and fipf — fipi = fi p . Combining these boundary condi- 
tions with the Landauer-Biittiker formula, the voltages 
V p (V p = Hp/e) in four voltage probes can be obtained, 
then the longitudinal resistance i?i4,23 = (V2 — V3)/Iu 
and Hall resistance i?i4,26 = (V2 — Vq)/Ii4 are calculated, 
here I14 = —I\ e = he- The resistances contain the prop- 
erties i?i4,26 = -^14,35 and i?i4 ; 23 = -^14, 65- 

Fig. 2a and 2b show the longitudinal and Hall resis- 
tances, iii.4,23 and i?i4 i 26 J versus the energy eo and the 
exchange split M at an external magnetic field <f> = 0.005. 
Due to the QSHE and QHE, both #14,23 and Ru^e may 
be non-zero, and they both exhibit plateau structures. 
The plateau values are determined by the filling factors 
and v^. For the fixed filling factors and t^, #14,23 and 
#14.26 maintain their plateau values regardless of eo and 
M. By considering the carriers transport along the edge 
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FIG. 2: The panels (a) and (b) show the resistances -Ri4,23 
and i?i4,26 (in the unit of h/e 2 ) vs the exchange split M and 
energy e for TV = 80 and (j> = 0.005. 
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states, the plateau values can be analytically derived;^ 3 - 
#14,23 = and #14,26 = [l/(^t + n)]V e2 for ( v tj v i) = 
(+, +) or (-,-), and #14,23 = [knl/(KI 3 
and #14,26 = sign{^)[(\^\ 2 - H 2 )/(KI 3 + 
for (y^,v^) = (+,— ) or (—,+). Some plateau values for 
low z/f , vi have been labeled in Fig. 2. The numerical re- 
sults in Fig. 2 are in excellent agreement with the analytic 
plateau values (the differences between them are less than 
1CP 6 ). Furthermore, #14,23 and #14,26 have the follow- 
ing properties: While | e 1 > \M\ with (z/f, z/jj = (+, +) or 
(—,—), the longitudinal resistance #14,23 is zero and only 
the Hall resistance #14.26 exists because the spin-up and 
spin-down carriers are simultaneously either electron-like 
or hole-like and move in the same direction. On the other 
hand, while |eo| < \M\ with (v^,v±) = (+, — ) or (—,+), 
the Fermi energy Ep is located between eq + M and 
£0 — M, the longitudinal resistance #14,23 emerges since 
now the spin-up and spin-down carriers move in opposite 
directions for a given edge, (i) While v-\ = —v\. = v, 
the Hall resistance #14,26 = 0, only the longitudinal re- 
sistance #14,23 exists with the value (l/2v)h/e 2 . This 
means that only the QSHE emerges and the QHE van- 
ishes in this region. In this case, the system has the CT 
invariance. Furthermore, if z/f = — = ±1, #14.26 = 
and #14,23 = (l/2)(h/e 2 ). Now the observed phenom- 
ena are completely the same with the QSHE from the 
SOIj£~— i 10 ] 11 but their physical mechanisms are different, 
(ii) While 7^ —v^ but still with (z/j-, v^) = (+, — ) 
or (—,+), #14,26 is now non-zero since the numbers of 
the spin-up and spin-down edge states arc different. In 
this case, both resistances #14,26 and #14,23 have non- 




FIG. 3: (Color online) The resistances -Ri.4,23 (a) and #14,26 
(b) vs the energy e for N = 80 and <j> = 0.005. 



zero quantized plateaus and the QSHE and QHE coexist. 
Fig. 3 shows the resistances #14,23 and #14,26 versus the 
energy eo for a fixed M (i.e. along the horizontal lines in 
Fig. 2), and it clearly shows that the quantum plateaus 
persist very well. 

Up to now, we demonstrate the existence of QSHE in 
the ferromagnetic graphene from both physical picture 
and detailed numerical calculations. In the following, 
we study the properties of the spin Hall resistance # s , 
a measurable quantity robust to dephasing 2 — and well 
reflecting the topological invariance of the system. We 
again consider the four-terminal Hall bar. But now the 
transversal terminals 2 and 4 are spin-biased probes with 
boundary conditions I p f = I p ± = (p = 2,4). Here the 
spin Hall resistance R s is defined as the transversal spin 
bias over the longitudinal charge current: R s = (p2-f — 
A 1 2i)/e/i 3 = ~(pit~ ^ii)/eh 3 . Since the spin bias n„ t - 
/j, n l is experimentally measurable, so is the # s i 25 i 26 Fig. 4 
shows R s versus the energy eo for different ferromagnetic 
exchange split M and magnetic field <j>. For |eo| > \M\ 
with = (+,+) or (—,—), # s = 0. On the other 

hand, while |e | < \M\ with (z^,n) = (+;~) or (->+)> 
# s exists. # s exhibits the quantum plateaus, and its 
plateau values are at [l/(|z/f|+ |^|)]/?,/e 2 . For a small M 
(e.g. M = 0.02i or 0.05i in Fig. 4b) or a high magnetic 
field <f> (e.g. = 0.005 in Fig. 4a), (v^, z/4.) can only equal 
to (1,-1), so only the plateau of # s = h/2e 2 emerges. 
But for a large M or a small magnetic field <fi, 
may be (1,-3), (3,-1), (1,-5), (5,-1), etc, then the 
plateaus of # s = h/4e 2 , h/6e 2 , etc, are also possible. 

Finally, we examine the disorder effect on the spin Hall 
resistance # s . Here we assume that the disorder only ex- 
ists in the central region (see dotted box in top right 
inset of Fig. la). Due to the disorder, the on-site energy 
£0 — crM for each site i in the central region is changed to 
to + Wi — aM , where Wi is uniformly distributed in the 
range [— W/2, W/2] with the disorder strength W. Fig. 5a 
shows # s versus the energy eo at the different disorder 
strengths W and Fig. 5b shows # s versus the disorder 
strength W at different energies eo- The results show 
that the quantum plateaus of R s are very robust against 
the disorder because of the topological invariance of the 
system. The quantum plateau maintains its quantized 
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FIG. 4: (Color online) (a) shows R s vs e for M = 0.05 and 
(b) shows Rs vs eo for <f> = 0.005. The parameter N = 80. 
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FIG. 5: (Color online) R s vs eo (a) and R a vs the disorder 
strength W (b) with the parameters N = 40, (j> = 0.007, and 
M — 0.05. The curves in (a) and (b) are averaged over up to 
1000 and 8000 random configurations, respectively. 



value very well even when W reaches 2 (see Fig. 5a and 



5b). Since the plateau is so robust and stable, its value 
can be used as the standard for the spin Hall resistance. 
In addition, even for a very large disorder strength W 
(e.g. W = 5 or larger), the plateau value only slightly 
decreases while maintaining the plateau structure (see 
Fig. 5b). This is because although the disorder strongly 
weakens the spin bias /i2t ~ ^24., it also weakens the lon- 
gitudinal charge current J13, so the value of R s is af- 
fected less. This means that in the large disorder limits 
(W — > 00), although the QSHE is broken, the SHE still 
holds. 

In summary, we predict a new QSHE in the ferromag- 
netic graphene film. Unlike the QSHEs studied so far, 
the origin of this QSHE is not caused by the spin-orbit 
interaction. The results also show that the system can 
exhibit the QSHE, the QHE, and the coexistence of the 
QSHE and QHE, depending on the filling factors of the 
spin-up and spin-down carriers. Due to the QSHE and 
QHE, both the longitudinal and Hall resistances exhibit 
the plateau structures. The plateau values (in the unit 
of h/e 2 ) are at 1/2, 1/6, 3/28, for the longitudinal re- 
sistance and at ±1/2, ±1/4, ±1/6, ±2/7, for the Hall 
resistance. In addition, the spin Hall resistance has also 
investigated and found to be robust against the disorder. 
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